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We present a quantization condition for the spectrum of a system composed of three 
identical bosons in a finite volume with periodic boundary conditions. This condition gives 
a relation between the finite volume spectrum and infinite volume scattering amplitudes. 
The quantization condition presented is an integral equation that in general must be solved 
numerically. However, for systems with an attractive two-body force that supports a two- 
body bound-state, a diboson, and for energies below the diboson breakup, the quantization 
condition reduces to the well-known Liischer formula with exponential corrections in volume 
that scale with the diboson binding momentum. To accurately determine infinite volume 
phase shifts, it is necessary to extrapolate the phase shifts obtained from the Liischer formula 
for the boson-diboson system to the infinite volume limit. For energies above the breakup 
threshold, or for systems with no two-body bound-state (with only scattering states and 
resonances) the Liischer formula gets power-law volume corrections and consequently fails 
to describe the three-particle system. These corrections are nonperturbatively included in 
the quantization condition presented. 

I. INTRODUCTION 

Quantum chromodynamics (QCD) as the underlying theory of strong interactions is responsible for 
a wide range of phenomena in nuclear and particle physics. Non-perturbative numerical evaluations of 
the QCD path integral through lattice QCD (LQCD) have started to bridge the gap between nuclear 
physics and the underlying theory of QCD. With ever- increasing computational resources and ever- 
decreasing uncertainties in LQCD calculations of various quantities, it is reasonable to expect few-body 
nuclear reaction cross-sections to be reliably evaluated in upcoming years. This program will allow for 
the determination of nuclear few-body forces in a model-independent fashion which will consequently 
lead to a reduction of the systematic uncertainties of many-body calculations. Reaching that point, 
however, requires developing the formalism to relate the physical observables of interest to the LQCD 
correlation functions that are evaluated in a finite Euclidean spacetime. 

For two-body systems below inelastic thresholds, the energy spectrum in a finite cubic volume 
with periodic boundary conditions fully determines the infinite volume scattering phase shifts, up to 
exponential corrections in the volume. These corrections are negligible provided that interactions are 
fully localized inside the volume. This idea, which is due to Liischer [H [2], has been extensively 
implemented in obtaining scattering parameters and binding energies of meson-meson, meson-baryon 
and baryon-baryon systems (see Refs. [HHS] for some recent works in this direction). This approach 
has also motivated studies of several hadronic resonances via LQCD as reviewed in Ref. [10]; in 
particular the two-body coupled-channel generalization of the Liischer formalism |llH15j allows one to 
access physics beyond inelastic thresholds where new two-particle channels open up. Going beyond 
two particles in this finite volume (FV) formalism, however, has been a major obstacle in the field. 
Recently several efforts have been initiated to address this problem via different approaches [16ft23j : 
aiming to make a connection between the elements of the S-matrix and the spectrum of the three- 
particle system in a finite volume. These benchmark theoretical calculations give one confidence that 
the infinite volume elements of the S-matrix can in fact be determined from the discretized spectrum 
in the finite volume even above the inelastic threshold, and has motivated us to further explore this 
connection using a non-relativistic field theoretical approach. 

In this article, we derive the quantization condition (QC) for the spectrum of a system composed of 
three identical bosons in a finite volume with periodic boundary conditions. The quantization condition 
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gives a relation between the FV spectrum and infinite volume scattering amplitudes, Eqs. (20), (19). 
This quantization condition in general must be solved numerically, since the relation between the FV 
energy eigenvalues and three-particle scattering amplitudes is not algebraic. We pay close attention 
to systems with an attractive two-body force that allows for a two-body bound-state, a diboson, and 
energies below the diboson breakup. For these theories in this energy regime, the quantization condition 
reduces to the well-known Liischer formula with exponential corrections in volume with a length scale 
that is dictated by the inverse diboson binding momentum. In other words, the boson-diboson scattering 
phase shifts can be obtained from the three-particle spectrum using the following relation 

p e~"''i^ 

cot 6Bd = 47r Coo(qS) + ' (1) 

where = ^J | (mE* + is the momentum of the boson and diboson in the center of mass (CM) 
frame of boson-diboson system with being the relative momentum of the two bosons in the diboson 
in the CM frame of the diboson, m is the mass of the three identical particles, E* is the CM energy, 7^ 
is the binding momentum of the diboson in the infinite volume limit, dsd is the scattering phase shift of 
the boson-diboson system, L is the spatial extent of the cubic volume, and r] is an unknown coefficient 
that must be fitted when extrapolating results to the infinite volume. Given that the diboson is bound, 
q'o' < and — ?■ —7^ as the volume goes to infinity, cf^ is a kinematic function that is given by 

J (27r)3 



k 



*2 



where V denotes the principal value of the integral, and for non-relativistic particles k* = k — aP. P is 
the total momentum of the boson-diboson system, and a = ^^^^ |21| : so for a diboson that is twice 
as massive as the boson a = |. In addition to exponential corrections that are governed by the size 
of the bound-state wavefunction, there are other exponential volume corrections to the above Liischer 
relation that are arising from the the off-shell states of the 2-1-1 system. These corrections however are 
subleading compared to the exponential corrections denoted in Eq. ([T]), and will be discussed in Sec. 
mil in more details. 

In order to reliably use such analytical formula one must necessarily be in the regime where 7rfL is at 
least 4 so that the infinite volume phase shifts of the bound-state particle scattering can be obtained with 
a few percent uncertainty. This is an important distinction compared to the the two-body problem where 
the dominant finite volume corrections to the Liischer formula scale like ~ e~^^^ where m-,^ denotes 
the mass of the pion. So although volumes of the order of 6 fm or greater would reliably recover, for 
example, tttt scattering phase shifts at the physical pion mass, in order to accurately recover phase shifts 
for deuteron-neutron scattering one would naively need L > 17 fm. But presumably upon quantifying 
the coefficients of these exponentials, linear combinations of these exponential corrections can be formed 
for different boost momenta of the three-particle system so that to cancel out the leading exponential 
corrections to the above quantization condition, and therefore reduce the size of the volumes needed for 
a reliable determination of the phase shifts to L > 12 fmPJ Note that the next to leading exponential 
correction due to the size of the bound-state scale as e~^'^"''^^ / L. The quantization condition shown in 
Section [TT| demonstrates that for energies above the diboson breakup Eq. ([T]) gets power- law corrections 
associated with new possible states that can go on-shell and the quantization condition must be solved 
numerically. 

Studying few-body systems in a finite volume is not only pertinent for the determination of nuclear 
reactions cross sections, but will also have an impact in our understanding of excited states and res- 
onances of QCD. This has recently motivated Guo et al. [23] to explore the three-body problem in a 
finite volume. The formalism developed in this article, however, clearly illustrates that the quantization 
condition for energy eigenvalues of the three-particle system in a finite volume in the isobar approxima- 
tion [25H28] that is derived in Ref . [23] , fails to describe the physics of the three-particle processes these 
authors have been interested in, for several reasons that are summarized in subsection III A| Basically 



^ For a discussion of the improvement of the volume dependence of deuteron binding energy in LQCD calculations see 
Ref. [24]. 



3 



although the quantization condition of Ref. |23j takes a simple algebraic form, there are power-law 
corrections in volume for resonant two-body sub-systems (which has been the motivation of Ref. |23j 
for studying the three-body systems) that invalidates the use of their result. These corrections arise 
from the fact that the two-particle sub-system is not compact but rather an extended object that, that 
regardless of its resonant nature, can sample the boundaries of the volume. These power-law corrections 
are nonperturbatively included in the quantization condition presented, Eqs. (20). The quantization 
conditions of Ref. [23] also do not incorporate the possibility of new on-shell three-particle channels 
that should turn the quantization conditions to a coupled-channels form. 

In addition, we comment on another source of systematic uncertainty in previous finite volume three- 
body calculations [T6HT9] that has not been fully addressed. This systematic uncertainty is due to the 
reduced symmetry of the calculations in a finite cubic volume with periodic boundary conditions. In an 
S-wave 2-|-l scattering where the two-body sub-system is also in the S-wave, the nearest partial-wave 
mixing, which is familiar from the Liischer formalism, occurs with the D-wave two-particle state for 
the identical particles [29]. This is due to the fact that in the CM frame of the three-particle system, 
the two-particle sub-system is necessarily boosted. As a result, the spherically symmetric S-wave dimer 
field which has been proven to be well suited for studying the three-particle system in the infinite 
volume |3UH32j . does not incorporate this additional feature of two-particle wave-function in a finite 
cubic volume. In other words, the D-wave partial wave mixing has been unaccounted for in previous 
numerical calculations by Kreuzer et al. |16H19j , which is the leading systematic error of their results. 

The rest of this article is organized as follows. In Sec. |IlJ after introducing the dimer formalism, 
we review the connection between the finite volume spectrum of the boosted system and the two-body 
scattering phase shifts. The advantages as well as limitations of the dimer formalism for this particular 
calculation are discussed. By means of the dimer formalism, the energy spectrum of three identical 
bosons with two-body and three-body interactions in the finite volume is determined from the poles of 
the finite volume three-body correlation function. The features of this quantization condition for the 
energy eigenvalues of the system are discussed. In particular, it is demonstrated that this condition 
gives an integral equation that relates the spectrum and the elements of the S-matrix. In Sec. Ill, the 
Liischer formula for the elastic scattering of a bound-state and a particle is deduced from the energy 
quantization condition. The systematic corrections due to the size of the composite object as well as 
those due to the presence of off-shell states are identified. We devote subsection III A to comparing the 
result of Ref. [23] to the result presented in this article, and will comment on its validity in studying 



three-body systems in the isobar model. We conclude in Sec. IV, and briefly discuss the strategies that 
a practitioner needs to follow in order to obtain the scattering quantities of the three-particle system 
from a LQCD calculation of the spectrum. 



II. THREE PARTICLES IN A FINITE VOLUME: QUANTIZATION CONDITION 

In studying a three-body problem, it is customary to divide the system of three particles to a system 
of two particles interacting in a given partial-wave J^, and a third particle, called the spectator, which 
interacts with the two-body system with angular momentum Jud- In particular, the diagrammatic 
representation of few-body scattering amplitudes is greatly simplified when using an auxiliary S-wave 
dimer field \33\ |33] which non-perturbatively sums all 2 — )• 2 diagrams. 

Consider three bosons with mass m and total energy and momentum (E, P) in the lab frame. The 
total CM energy of the three-particle system, E* is then given hy E* = E — Also the relative 
momentum of the spectator boson and the dimer in the CM frame of three particles, q*, is related to 
the momentum of the spectator boson in the lab frame, q by q* = q — ^. The total CM energy of 
boson-dimer system can be written as E* = q*'^ /m + 3g*^/4m, where q* is the relative momentum of 
the two bosons inside the dimer. 

One can write the full dimer propagator in the infinite volume by summing up all the 2 — )• 2 
interactions as depicted in Fig. ([T| a). By requiring the infinite volume two-body scattering amplitude 
to be recovered from the full dimer propagator, one obtains 

rV^{E,,cf) = ^ , (3) 

q* cot Od — iq + «e 





V'^= = + 

(a) 

v''= = = + 

(b) 

FIG. 1: a) Diagrammatic equation satisfied by the full dimer propagator in a) infinite volume and b) 
finite volume. The grey (black) band represents the full infinite (finite) volume propagator while the 
double lines represent the bare propagator. 



where E2 = E* — q*"^ /2m is the total energy of the dimer. 5d denotes the S-wave scattering phase shift 
of the two-boson system, and r is its effective range. Similarly one can work out the finite volume case, 
Fig. ([T| b), where due to the periodic boundary condition, the lab frame momenta are discretized, 
q = where n is a vector with integer components, and the loop integrations are replaced by the 

corresponding sums, J ^^^^3 — )• -p^ Then it is straightforward to show that the full finite volume 
dimer propagator can be written as 

^^^(i?2,q*) = ""^'P . , (4) 

q* cot 5d - 47r clQ{q*^+ie) + ie 

where the kinematic function is defined in Eq. For the case of a dimer composed of identical 
bosons, one therefore has k* = k — 

The spectrum in a finite volume can be obtained from the poles of the two-particle propagator or 
equivalently from the poles of the finite volume dimer, Eq. Q, 

CcotJrf = 47r cgo(C^) ' (5) 

where is the At*'' solution to the quantization condition for a boosted two particle system [JHl |35] . 
As will be discussed in great details, these poles play an important role in the three-body sector and 
will be referred to as Liischer poles. Note that this result is equivalent to the non-relativistic limit of 
the result obtained in Refs. [29l [35l [36] for the boosted systems of particles with identical masses, and 
for that of systems with unequal masses [Ml [SZl [38] . 

The simplification due to the use of the dimer field comes at the cost of truncating the precision of 
the effective range expansion of the two-body sector at next to leading order (NLO). This, however, 
can be shown to be systematically improved by including higher order terms in the effective range 
expansion in constructing such a dimer field. The other systematic error associated with introducing 
such an auxiliary field in the finite volume arises from the fact that the spectrum that is obtained 
by looking at the poles of the dimer propagator corresponds to two bosons in an S-wave, and does 
not incorporate the partial-wave mixing. This partial-wave mixing is expected since the eigenstates of 
Hamiltonian in a finite cubic volume do not respect the full symmetries of the rotational group, and 
are instead identified by the irreducible representations (irreps) of the Cubic (octahedral) group. The 
irreps of the cubic group are in general reducible under the irreps of the full rotation group. So, for 
example, the ground-state of two bosons in a finite cubic volume is an eigenstate of the Ai irreducible 
representation of the cubic group, which in the CM reference frame has overlap with not only J = 
but also J = 4, 6, 8, . . . [H [2]- For two identical particles moving with a non-zero total momentum, 
the symmetry is further reduced and the ground state has overlap with J = 0, 2, 4, 6, 8, . . . [291 [35l [36]. 
Although, this mixing between higher partial-waves is explicitly manifested in the pole structure of 
the two-particle propagator [H [21 l29l l35l |36] ; it is not present in the dimer formalism. Nevertheless, 
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(b) 

FIG. 2: a) Diagrammatic expansion of the three-body correlation function in the finite volume. 
Kj, denotes the three-body Bethe-Salpeter kernel and A^, (^3) is the overlap the creation 
(annihilation) dimer-boson interpolating operator has with the initial (final) state with total energy E 
and total momentum P, b) The effective three-body Bether-Salpeter kernel, K^, is composed of a 
three-body contact interaction as well as two-body contact interactions via the exchange of a single 
boson. 



this technique has proven to be advantageous for studying three-body physics in both infinite and 
finite volumes [T6HT91 [311 [321 1391431] , and as a first step toward solving the three-body problem in the 
finite volume, it is reasonable to start with this dimer field. One needs however to keep in mind the 
systematics of the calculation as presented. 

In order to determine the energy eigenvalues of the three-particle system in a finite volume, one can 
solve for the poles of the three-particle correlation function as depicted in Fig. ([2| a) . Algebraically, 

CY{E,V) = E ^3 (qi) ^V''{E - - Qil) 

<=ii 



1+ E n ( ^E^^3(qn-i,q„;P,i?)^^^(i?- £,|P-q„|) ) A'3 (q„0 



n.'=2n=2 \ q„ 



, (6) 



where A'^ (^3) is the overlap the annihilation (creation) dimer-boson interpolating operator, as (cj), 
has with the initial (final) state with total energy E and total momentum P. Note that we have 
suppressed the total energy and momentum dependence of the overlap factors in our notation. The 
interactions between three bosons are incorporated in an effective three-body Bethe-Salpeter Kernel, 
Ks, Fig. ([2} b), 

iKsip, k; P, E) = -^9s - —^-^^^^z^z^— , (7) 

^ 2m 2m 2m + 

where the incoming (outgoing) boson has momentum p (k) and the incoming (outgoing) dimer has 
momentum P — p (P — k), and {E,F) denote the total energy and momentum of the three-particle 
system as before. The first term in the Kernel, Eq. ([7|), is the three-body contact interaction, while 
the second term describes the interaction of three particles via exchange of an intermediate particle 
through two-body contact interactions. 

The finite volume contribution to the first term in the expansion of the three-body correlation 
function, Eq. ([g]), can be evaluated easily using the Poisson resummation formula and the kinematic 
relations between the CM and lab frame momenta as presented earlier. 



zsE-/ |^U3(qi)^2?^(i?-|^jp-qii)4(qi) = E^3(9:)i5aT(5M(g:)- («) 

Ql I K 



2 ^'^ 




= ^^^^^^^ + 

FIG. 3: The NLO contribution to the three-body correlation function with one insertion of the 
three-body kernel. 




While ^3 and A'^ in the LHS of Eq. (jSj) are functions of the relative coordinate qi , they are represented 
as vectors in the space of the boson-dimer angular momentum, Jb^, in the RHS, and are evaluated at 
the poles of the dimer propagator, q*, and the sum runs over these poles. 60^ is a matrix in the same 
angular momentum basis whose elements are defined by 



imi,/2»Ti2 



(9) 



with 



imi,/2ni2 



. "t q^ 



+ i 



Im 



'Att 



(10) 



The kinematic function cf is defined in Eq. with a = i, since the dimer is twice as massive as the 



boson. The on-shell CM momentum of the boson-dimer system, g*, is defined by k*'' pole of the FV 
dimer propagator, g*^ = mE* — and i?^ is its residue at the pole. Explicitly, 



lim iV'^iE 



2m' 



q| 



iR 



V 



iRl 



3q 



*2 



(11) 



where 



R 



V 



mr 



d 



dq 



*2 



q* cot Sd — An q 



|P-q| 

00 



-*2 



=9k 



(12) 



Note that the poles of the FV dimer propagator correspond to the energy eigenvalues of the boosted 
two-particle system in the finite volume, Eq. i.e. the Liischer poles. 

Equation Q reflects the fact that, unlike the two-body case, a single on-shell condition does not 
simultaneously fix the relative momenta of the boson-dimer pair as well as that of the bosons inside 
the dimer. What is important to observe here is that for a given CM energy E* , there is only a finite 
number of "channels" Ne* that can go on-shell, each being identified by a particular configuration of 
the boson-dimer relative momentum and the relative momentum of the two bosons (of the dimer), 




lN,,,\i^{mE*-q*j^^J 1 ;>(13) 



These channels contribute to the quantization condition since q'^/m < E*. This observation makes 
the analogy to the "coupled-channel" systems self-evident. For CM energies that are below the dimer 
energy, E* < ql^/m, the energy is not sufficient to allow the three-particle system to go on-shell. 
Subsequently, these states can be neglected as they give rise to exponential corrections in volume 
rather than power-law. Furthermore, similar to the two-body case, the on-shell condition does not 
fix the directional degrees of freedom of the relative momentum of the 2-1-1 system, and therefore 
it is convenient to upgrade all finite volume quantities into infinite-dimensional matrices in angular 
momentum. 
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The calculation of the second term in the expansion of the correlation function, Eq. Q, is more 
involved as it comes with one insertion of the three-body Kernel, Fig. ([s]), and due to the one boson 
exchange contribution couples the momenta running into the loops. 



qi.q2 



-m 



E 



qii 



2m 



2m 



2m 



+ 26 



2m 



|P-q2lK3(q2). (14) 



Although at the first glance, there appears to be poles arising from the exchange boson propagator, 
one can verify that the poles of the three-body kernel are exactly canceled by the zeros of the full finite 
volume dimer propagatoi[^ As a result, the only power law volume dependence of such diagrams arise 
from the poles of the dimer propagator only. Given this observation, it is straightforward to show that 



CUE) 



2m 



c^^\)Y,[KMl,<ll;E'')6gl{ql)AM 



Ne* 

E 

k,k' 



|P-q2lK(q2) 



(15) 



where a summation over angular momentum is understood for the terms inside the brackets. The 



summation over the two-body Liischer poles is left explicit. The result of Eq. (15), along with the fact 
that the dimer propagator can be decomposed in a series over its poles, 



W^{mE -3q'^/A,q) = ^ 



iRl 



7*2 



+ ie 



(16) 



suggests that the dimer propagator can be upgraded unto a diagonal matrix in the space of A''^;* available 
FV states which is a useful representation when performing the sum over all diagrams contributing to 
the correlation function. Each element of this matrix is then effectively a single particle propagator 
with the corresponding FV pole and residue that contain finite volume dependence of the propagators, 



[iV^{mE -3q^/A,q)]^ 



iRl 



q*2_ 



1*2 



+ ie 



(17) 



Given the simplifying feature of the FV loop sums as is evident from Eq. ( 15 ), and the representation 



of the FV dimer a matrix over available channels, Eq. ( 17), it is straightforward to sum over the infinite 
number of terms appearing in the boson-dimer correlation function, Eq. ([6]). Denoting the boson-FV 
dimer propagator as Qy, one can show that Eq. ^ is equal to 



C^{E,P)-C^y{E,P) = iA3 



(1 - g^M^)6g^ 



1 



^{i-M^g^) 



A' 



3 ) 



(18) 



where C^y{E,P 



My" is defined as the sum over all infinite volume di- 



iA:,g^{l-M^g^)A',. 

agrams containing a boson and a finite volume dimer, Fig. (|4| c), and can be interpreted as the 
non-renormalized infinite volume scattering amplitude between a boson and a FV dimer]^ That is to 



^ This important observation was first pointed out to us by Michael Doring and Akaki Rusetsky for the relativistic 
three-particie systems |45| . 

^ The difference between renormalized and non-renormahzed scattering ampiitudes wili be expiained shortiy. 
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(a) 




(b) 




(c) 



FIG. 4: a) Diagrammatic representation of the inhomogenous integral equation satisfied by the 
three-body scattering amphtude in the infinite volume, b) The corresponding sum equation satisfied 
by the FV scattering amplitude, c) The diagrammatic representation of the integral equation satisfied 
by7W~, Eq. 



say, while the relative momentum between the dimer and boson is continuous in Aiy, the relative 
momentum of the two bosons inside the dimer remains discretized, 

j^M^P, q; P, E)5V''iE - 1^, |P - q|)>^y (q, k; P, E) 

(19) 

where 6V^ =V'^ -V°° and is the non-renormalized infinite volume scattering amplitude, Fig. 

(j4| a). For comparison, the full FV scattering amplitude of the dimer-boson system is also depicted in 
Fig. ([4| b), where all relative momenta between the three-particles are necessarily discretized]^ Note 
that the boson-dimer propagator and ^A'y are evaluated on-shell where the scattering energies of each 
boson-dimer channel is given by (f^ jm. 



The poles of the FV correlation function, Eq. (18), determine the spectrum, 
Det(l + >(??(^g^) = detc 



detpw(l+A^^<5a^)J =0 , (20) 

where the determinant detoc is over the A'^^;* open channels and the determinant detpw is over the 
boson-dimer relative angular momentum. In practice it is necessary to perform a truncation over 
the partial-waves and choose a maximal angular momentum. This quantization condition however 
incorporate the partial- wave mixing due to the reduced symmetry of the boson-dimer wavefunction in 
the finite cubic volume as will be discussed in more details in the next section. 

The reason that the scattering amplitude quantities introduced above are not renormalized is that 
unlike single particle operators in a non-relativistic field theory, the dimer field corresponds to an 
interpolating operator that has overlap with two-particle states, and as a result must be renormalized 



* One should however note that due to the absence of asymptotic states in a finite volume, the interpretation of this 
quantity as the scattering amplitude is ambiguous, and it is only introduced as a useful mathematical quantity in 
analogy with the infinite volume scattering amplitude. 
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[401 H6] . The renormalization factor in the finite volume can be obtained from the residue of the FV 
dimer propagator, 



d 



dE* 



iV^ E 



2m' 



m 

w 



(21) 



.Ml 

4m 



Upon renormalizing the dimer field, tlierefore, one arrives at the normalized scattering amplitudes in 
the finite volume, e.g. (M^)^^' = {ZX)^/\M^)^^>{ZX,)^/\ 



The quantization condition, Eq. (20), resembles that of the two-body coupled- channel systems as 
presented in Refs. [21 [TB]. As discussed earlier this illustrates that a single on-shell condition does 
not fix the magnitude of both relative momenta and there is a freedom in scattering in any of finite 
number of available channels. The other characteristic of Eq. (20) is that it does not still provide 



a algebraic relation between the infinite volume scattering amplitude and the energy eigenvalues of 
the boson-diboson system, simply because A4y still has possibly large volume corrections arising from 
FV dimer propagator, Eq. (19). Despite all these complexities, this quantization condition not only 



gives a better insight into the three-body problem in a finite volume, it automatically reduces to the 
Liischer quantization condition for the bound-state-particle scattering below the bound-state breakup, 
up to exponential corrections that are due to the size of the bound-state wave-function. This will be 
discussed in the next section in more details. 



III. BOSON-DIBOSON SCATTERING BELOW THE BREAKUP THRESHOLD 

The formalism developed in the previous section does not assume any specific form for the interac- 
tions in the three-body system. Therefore the result presented is universal regardless of the nature of 
the interactions or whether the theory contains any number of two-body or three-body bound-states. 
In this section, though, we consider a theory with an attractive two-body force which allows a two-body 



bound-state, a diboson. We will show how Eq. (20) reduces to the well-known two-body result below 
the diboson breakup. 

For such energies there is only one state that can go on-shell and introduce power-law volume 
corrections, the boson-diboson state. By restricting to /max = 0, the low-energy parametrization of the 
scattering amplitude becomes that of a two-particle system in an S-wave with masses mi = ^ = m, 

Stt 

m cot Ssd - iqo 



^Bd = — ^,_,,__ (22) 



where q^^ /m = |(i?* — q'^ /m) is the boson-diboson scattering energy in the CM frame, q^ jra is 
the boosted diboson FV binding energy, and bsd denotes the boson-diboson scattering phase-shift. 



However, this is the non-renormalized quantity My that appears in the QC, Eq. (20), and not the 
physical scattering amplitude. Here we argue that by introducing a systematic error of the order of 
eT'^'i^ jh to the final result, Eq. (Ilj), the scattering phase shifts can be derived from the QC, Eq. (20) 
after replacing bQ^ with bQ^M^ = SG^Msd- Id denotes the infinite volume binding momentum 
of the diboson which satisfies, 

(rcot5,-ir)l^*=,^, =0. (23) 

The first step to prove this claim is to note that the bound-state pole of the FV dimer propagator is 
exponentially close to the bound-state pole of the infinite volume dimer, = + 0{e~"''^^ / L), which 
is evident from Eq. ([s]) after analytically continuing the momentum to the imaginary axis. Now in 
evaluating My , one needs to perform a series of coupled integrals as is given in the first line of Eq. 



(19). For negatives energies, the only singularity of the integrands in the range of integration occurs 
when the diboson pole of the FV dimer propagators is reached. The contribution to the integrals due 
to this singularity is proportional to the residue of the FV dimer at that pole. Since the residue of the 
infinite volume dimer propagator at the diboson pole, 

mr 



T?00 

Rd 2 



d 



dq 



^2 {q* coi5d-iq*)\q*2^_^2 



(24) 
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is exponentially close to its FV counterpart, Eq. pilpj 



Rf [1 + 0(e-^^^) 



"^Ve jj , (25) 

one can replace with 1)°° up to the exponential corrections that scale by the size of the bound-state 
wave-function. Consequently, from Eq. (19) one observes that M.^ is equal to Aiy up to exponentially 
small corrections. Note that A^?? and are renormalized differently, however, the finite volume 



dimer field renormalization factor , Eq. (21) is exponentially close to the renormalization factor of 



the infinite volume diboson field Z°° around the bound-state pole, which is defined as 



d 



dE* 



iV°^ E 



r_ IP 

2m' 1-^ 



q| 



m 

TJOO 



(26) 



60^ in Eq. 



(20), with 50^ M'^ for elastic srocesses 



Therefore one can approximate 5Q^ M'^ 
that occur in this energy regime up to these exponential corrections as stated. 

Keeping in mind these exponential corrections, one can now apply the expression for the scattering 
amplitude in Eq. ([22]), to Eq. ([20|. Using Eqs. ([lO]), (|9j), one can recover the two-particle quantization 
condition for two-particle systems up to the exponential corrections explained above as shown in Eq. 
([1]) and reiterated here for clarity 

p-ldL 



ql cot 6Bd = 47r Coo(gS) + r/- 



L 



(27) 



This result confirms the postulation and numerical verification made by Bour et al. |22j that upon 
subtracting off the FV binding energy of the bound-state from the total energy of the three-particle 
system, the scattering energy eigenvalues of the bound-state-particle system can be reliably related to 
the scattering phase shift of the system through the use of Liischer formula for two-body systems after 
extrapolating to the infinite volume limit. This offers the practitioner a reliable method to extract the 
infinite volume phase shift of elastic bound-state-particle scattering by fitting to an exponential form. 

This result also illustrates that in order to obtain the boson-diboson scattering phase shift, not 
only one needs to determine the boosted three-particle energy spectrum, but also needs to obtain the 
scattering parameters of the boosted two-particle system first. It is also evident that if the interactions 
support a boson-diboson bound-state, a triboson, after analytically continuing the momentum in Eq. (fll) 

can be 



4m 



to the imaginary plane, Qq = i^Bdi the binding energy of the three-particle system, B3 
obtained easily via Eq. ([T]), as is well-known for bound-states appearing in the two-body sector [HEIIIT]. 
Alternatively, one can also solve for the triboson poles of the FV scattering amplitude from the FV 
counterpart of the STM equation. Fig. (|4] b), as is pursued in Refs. jT6Hl9]. 

The boson-diboson QC, Eq. ([T]), is a low-energy approximation of Eq. (20), which at NLO has two 



sources of exponential corrections. First, the QC receive corrections associated with the finite volume 
binding momentum of the diboson which scale like 0{e~^'^'^^ / L) at next to leading order. It also 
acquires exponential corrections associated with the truncation of off-shell states appearing in the de- 



composition of the dimer propagator, Eq. (16), as mentioned before. More explicitly, the next excited 
state of the boson-diboson system corresponds to a CM scattering energy of /m = |(£'* — qf^/m), 
where q^ is the boosted momentum for an unbound two-boson system. For E* < g^^/m, the three- 
boson scattering energy, ql'^/m, is negative which leads to exponential corrections of O (e~''i^/L) to 
the single-channel QC, Eq. M, which however are subleading compared to O (e~''o^/L) ~ 0{e~^'^^ / L) 
corrections. For sufficiently high energies, these exponential corrections become power-law in the vol- 
ume, and one necessarily has to study a coupled-channel system made up of a boson-diboson state and 



a three-boson state. For energies just above the diboson breakup Eq. (20) can be written as 



1 + 7W°^ 



V,Bd- 



'Bd ^^Bd 



1 + M°^ 



V,BBB~BBB 



5Q 



BBB 



\M 



00 

V,Bd- 



-BBbI"^ ^^Bd ^GbBB 1 



(28) 



There is another correction to the residue function at the diboson pole that occurs at O [e~'^''^ /jdL) . Since for 7dL > 1 
the diboson does not fit in the volume, and the finite volume formalism is no longer valid, one must make sure to use 
sufficiently large volumes for shallow bound-states so that -jdL ^ f . ft then follows that these corrections are subleading 
compared to the 0{e corrections in Eq. (25 1 and could be ignored. 
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where SQ^^^ and SQl^^^ are respectively the boson-diboson and three boson propagators, A^y^_^/ 

denotes the elements oi Aiy for the k^^ (k'*^) initial (final) state. For such energies, the approximations 
made before are no longer valid and determination of infinite volume scattering cross sections from the 
finite volume spectrum requires numerically solving an integral equation for A4y, Eq. (19). 

Lastly we comment on the systematic uncertainties associated with the dimer formalism and partial- 
wave mixing. Assume, for example, that both the dimer and the boson-dimer wavefunctions are pro- 
jected onto the ^4^*" irreducible representation of the cubic group and that the three particles are degen- 
erate. Then in the boson-dimer CM frame, the system has an overlap with ( J^, JBd) = (0, 0) as well as 
( Jrf, JBd) = {(2, 0), (4, 0), (0, 4), (2, 4), (2, 6), . . .} angular momentum states, with the leading contami- 
nation arising from the D-wave dimer. As discussed in Sec. |ITJ this is due to the the fact that the dimer 
in this 2 + 1 body set-up is boosted and its symmetry group in its CM frame is reduced compared to 
the original cubic group [29] . If one then proceeds to consider a reference frame where the dimer-boson 
system has non-zero momentum, then the ground state will have overlap with (Jrf, Jed) = (0,0) as 
well as ( Jrf, Jsd) = {(0, 1), (2, 0), (2, 1), (0, 2), (2, 2), . . .} angular momentum states. This is because the 
dimer-boson is effectively a two-particle system where one of the particles is twice as massive as the 
other, and therefore S and P-wave mixing is unavoidable . As a result one needs to simultaneously 
determine S and P-wave scattering parameters. Note that although the dimer field used in this paper 
is an S-wave filed which does not lead to inclusion of higher partial-waves in the two-body QC, the 
boson-diboson scattering QC, Eq. ([T]) fully incorporates the partial-wave mixing in the space of the 
boson-diboson angular momentum states. 



A. On the quantization conditions derived by Guo, et al. in the isobar approximation 

Developing a finite volume multi-particle formalism not only is crucial in studying nuclear reactions 
from LQCD, but also is a necessary step towards reliable LQCD studies of higher excited states as well 
as hadronic resonances. With this latter goal in mind, authors of Ref. [23] have attempted recently 
to present such a FV formalism for three-body hadronic processes that have a dominant two-particle 
resonant channel. Such a study is motivated since, as an example, the decay of J /ip to three pions has 
been shown to be well- approximated by decaying to a pion as well as the p resonance where the p can 
subsequently decay into two pions |l8] . Modeling a process that involves a three-particle final state via 
such a quasi two-particle intermediate state is called the isobar approximation [25H28j . and has been 
used by Guo et al. to study the spectrum of the three-particle system in the finite volume. Here we 
argue that the result presented by these authors fails to describe the three-particle processes that can 
be studied via the isobar model in the infinite volume, and by comparing their quantization condition 
with that of presented in this paper, confirm the validity of their result in one particular case that is 
different from processes which the isobar model concerns. 

Before discussing the result of Ref. |23], it is necessary to briefly review some important features 
of a resonance that makes it challenging to be extracted via a LQCD calculation. As is well known, 
a resonance is not an isolated eigenstate of the hamiltonian but rather a resonance in multi-particle 
scattering channels that have the same quantum numbers as the resonance state. In the limiting case 
where the resonance has an infinitely narrow width, one however recovers an isolated state of the 
theory. The key is that the resonance itself should not be confused with a composite object, e.g. the 
p meson can not be interpreted as a two-pion state; as soon as the dynamical coupling of p to the 
two-pion state is turned off, it no longer can decay to two pions. That is a key difference between a 
resonance that dominates the two-particle scattering amplitude, and a two-particle bound-state. As a 
result, the extraction of two-body bound-states follows directly from the identification of the negative 
interaction energy levels of the two-particle spectrum, and its binding momentum in the finite volume 
is directly related to the scattering phase shift via the standard two-body Liischer formula. To extract a 
resonance state from a LQCD calculation, however, requires mapping out the scattering phase shifts of 
the corresponding two-particle scattering channel as a function of energy using the Liischer formalism 
(for a recent implementation of this technique see Ref. ^49]). In addition, since most of resonances 
lie above the inelastic threshold, their extractions require studying a coupled-channel multi-particle 
scattering system. 
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With this difference in mind, let us discuss the result derived by Guo, et al. for three-particle systems 
in a finite volume within the isobar approximation. Based on a Hamiltonian formalism approach, Guo 
et al. presented two sets of Liischer formula that are derived from a relativistic Lippmann-Schwinger 
equation. One of these formulae describes the scattering of the two particles that compose the isobar 
in terms of the invariant mass of the isobar in the 2 + 1 system. This equation corresponds to Eq. ([s]) 
of this article where the scattering phase shift of the two-particle sub-system is related to the relative 
momentum of the particles in this sub-system. The other equation however, is effectively another two- 
particle Liischer formula describing the scattering of an isobar off the spectator particle. This relation 
offers an algebraic relation between the phase shift of the particle-isobar scattering system in terms of 
the relative momentum of the isobar and particle in the finite volume. Note that the relative momentum 
of the isobar and particle is correlated with the solutions to the Liischer equation for the isobar system 
at any given total energy which is a direct consequence of three-particle kinematics. 

There are two features of these quantization conditions that make them unrealistic solutions to the 
three-body problem within the model considered. First of all, as is discussed in detail through this 
article, for any given CM energy of the three-particle system, there are more than one configuration 
of momenta in the finite volume, see Eq. (13), that can go on-shell, and the corresponding channels 



need to be included in the quantization condition simultaneously. Although the authors of Ref. |23j 
point out that there are multiple combinations of relative momenta for each given energy that could 
satisfy their quantization conditions, their results fail to incorporate this coupled-channel feature of the 
problem. As a direct consequence, not only systematic errors due to nearby off-shell states of the 2 + 1 
system can not be addressed by their derivation, but most importantly their quantization conditions 
fail to describe the true behavior of the system as soon as an extra on-shell state is available at a given 
energy. It is important to reiterate that these "states" are nothing more than the kinematically allowed 
configurations of the three-particle system. 

The other misleading feature of these quantization conditions is that although the formalism is 
presented in the framework of the isobar model, no connection to the resonant nature of the two- 
particle sub-system has been made. Most importantly, the isobar, although introduced as a sharp 
resonance, has been treated as a composite system of two-particles. This composite system then has 
been studied via a two-body Liischer formula in a finite volume similar to two-body bound-states. The 
question that arises is that if the resonance is being treated as an isolated particle that scatters off 
another particle in this setup, which is a reasonable assumption in the limit of infinitesimal width, 
how could its energy spectrum be related to the scattering parameters of the two-body sub-system? 
On the other hand, as soon as a dynamical coupling to the two-particle state is turned on (which is 
the case for systems with three-body final states), the resonance is no longer an isolated state of the 
hamiltonian and its energy is not given by the standard Liischer formula, in contrary to the two-body 
bound and scattering states. One needs to take into account this distinction for one primary reason. A 
resonance state can get large volume corrections in a finite volume, and as a result, the correspondence 
between the isobar-particle spectrum in the finite volume and the isobar-particle scattering amplitude 
through a simple Liischer formula is ambiguous. This is analogous to our discussion of the 2 + 1 systems 
above the breakup where the FV two-particle state receives large power-law volume corrections and 
the boson-FV dimer scattering amplitude can no longer be approximated by the three-particle infinite 
volume scattering amplitude. That being said, our formalism demonstrates how one can recover the 
infinite volume scattering amplitude from the spectrum, but by no means should one expect an algebraic 
Liischer quantization condition for such a quasi 2 + 1 system. 

There is however one limiting case where the result of Ref. [53] is valid. As was discussed in 
Sec. in, the scattering phase shifts of boson-diboson scattering can be extracted from an algebraic 
Liischer quantization condition, up to exponential corrections with the leading piece scale by e~'^'^^ /L, 
Eq. This is a valid description of the process since in the energy regime where the two-body 

bound-state can not break up, there is only one on-shell momentum configuration that could lead to 
power-law corrections in volume. The quantization conditions presented by Guo et al. is therefore 
valid for theories that allow a deeply two-body bound-state, in order for the exponential corrections to 
be negligible, and in a energy regime where the bound-state can not disintegrate into its constituent 
particles. Then effectively one has an elastic scattering between two particles that can be studies via 
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a single- channel Liischer quantization condition, where the energy of one of the particles (the bound- 
state) is determined from another Liischer formula. It is only for these particular conditions that the 
equations given in Ref. f23j are valid. So in conclusion, although the quantization conditions of Ref. 
[23] could be applicable to deeply bound-state particle scattering below the breakup threshold, it cannot 
make a valid connection between the scattering parameters of an isobar model in infinite volume and 
the energy eigenvalues of three particles in a finite volume. 



IV. SUMMARY AND CONCLUSION 

Determining nuclear reaction cross sections directly from the underlying theory of QCD will impact 
our understanding of a wide range of phenomena. Currently LQCD is the only reliable way to carry out 
such ambitious program provided that the formalism that connects the physical scattering amplitudes to 
the lattice multi-nucleon correlation functions is put in place. With this in mind, we have determined a 
model-independent representation of the quantization condition for energy eigenvalues of three identical 
bosons in a finite volume with the periodic boundary conditions. Using a non-relativistic EFT, the FV 
three-particle spectrum has been shown to be related to the infinite volume S-matrix elements. For 
arbitrary energies, this correspondence requires solving an integral equation. With nuclear systems in 
mind, close attention is paid to scalar theories that support a two-body bound-state. It is shown that 
for energies below the diboson breakup, the quantization condition reduces to the well known Liischer 
result for two particles with unequal masses, with exponential corrections dictated by the size of the 
diboson, Eq. ([T]). Although physically intuitive, this is a non-trivial observation that will require the 
Lattice practitioner to extrapolate the phase shifts obtained from the Liischer formula for the particle 
bound-state elastic scattering to the infinite volume. 

In deriving the quantization condition for the three-boson system, an auxiliary dimer field can be 
used for convenience. This dimer field has the same quantum numbers as of the two-boson system in an 
S-wave, and enables one to study the three-body system as an effective 2+1 system. This simplification 
comes at the cost of neglecting the partial-wave mixing associated with a cubic finite volume with the 
periodic boundary conditions. By evaluating poles of the full three-boson correlation function in the 
finite volume, the quantization condition for the three-particle system with non-zero total momentum 



is derived. The poles are given by a determinant condition, Eq. (20), where the determinant is taken 
over dimer-boson relative angular momentum states as well as Ne* available boson-dimer eigenstates 
for each CM energy, E* . As is shown, the corresponding quantization condition has strong parallels 
with two-body coupled-channel systems |14^ 115] . However one has to be careful that this is not the 
physical scattering amplitude that directly shows up in the QC, but rather the scattering amplitude of 
a boson-FV dimer system. These two quantities are related to each other through an integral equation. 



Eq. (19). This manifestly shows that the determination of the scattering quantities of the three- 
particle system from the energy eigenvalues does not simply follow from two coupled Liischer formulae 
for two-particle sub-system, as well as 2+1 scattering system, as suggested in Ref. [23J. It has been 
demonstrated that as the energy of the system increases, more three-particle states can go on-shell, 
and as a result large FV correction due to the size of the two-particle sub-system can not be negleted. 



This feature is clearly embedded in the QC presented, Eqs. (20), (19). 

Furthermore, as is explained in detail, the exponential volume corrections from the off-shell excited 
states of the dimer are accounted for in the full quantization condition. For sufficiently high energies, 
these exponential corrections become power-law in volume and can no longer be neglected. Then one 
would have to consider a coupled-channel system where the number of channels are determined by 



the total CM energy of the three-particle system, as shown in Eq. (28) for energies just above the 
diboson breakup. The formalism presented considers three-particle with non-zero total momentum 
which eventually allows for more independent measurements at a given energy. This however leads to 
a practical complexity as the symmetry of the system is reduced, and the ground state of the system 
is expected to mix with the P-wave scattering state |2H [Mj [37] . The quantization condition derived 
predicts this mixing between S and P partial- waves, and indicates that the truncation of the determinant 
condition at S-wave could, in practice, introduce large systematics to the calculation. 

With these observations at hand, we argue that future LQCD studies of nuclear reactions and 
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resonances involving three-particle states will require the following steps. First, one needs to reliably 
determine scattering phase shifts for the two-particle sector from which one can obtain the boosted 
Liischer poles as a function of the boost momenta and energy. From there, one would proceed to 
obtained the three-particle spectrum. This requires high statistics to obtain multiple states with clean 
signals. Also in order to disentangle the coupled-channel nature of the three particle system, these 
calculations need to be performed with different boosts and in different volumes. In addition one has 
to simultaneously determine energy eigenvalues of three-particle states in different irreps of the cubic 
group to correctly deal with the partial wave mixing which is more severe than the two-particle sector. 
All of thes information should be simultaneously fitted to numerically solve the quantization condition 
presented. This would lead to an accurate determination of the three-body Bethe-Salpeter Kernel, 
which encodes all of the infinite volume physics up to the four-particle inelastic threshold. 
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